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We show that, in contrast to known results in the massive case, a minimally gauged 
massless Rarita-Schwinger field yields consistent classical and quantum theories, with a gen¬ 
eralized fermionic gauge invariance. To simplify the algebra, we study a two-component left 
chiral reduction of the massless theory. We formulate the classical theory in both Lagrangian 
and Hamiltonian form for a general non-Abelian gauging, and analyze the constraints and 
the Rarita-Schwinger gauge invariance of the action. An explicit wave front calculation for 
Abelian gauge fields shows that wave-like modes do not propagate with superluminal veloc¬ 
ities. The quantized case is studied in gauge covariant radiation gauge and T 0 = 0 gauge 
for the Rarita-Schwinger field, by both functional integral and Dirac bracket methods. In 
'ho = 0 gauge, the constraints have the form needed to apply the Faddeev-Popov method for 
deriving a functional integral. The Dirac bracket approach in To = 0 gauge yields consistent 
Hamilton equations of motion, and in covariant radiation gauge leads to anticommutation 
relations with the correct positivity properties. We discuss relativistic covariance of the 
anticommutation relations, and of Rarita-Schwinger scattering from an Abelian potential. 
We note that fermionic gauge transformations are a canonical transformation, but further 
details of the transformation between different fermionic gauges are left as an open problem. 


I. INTRODUCTION 


A. Motivations and Background 


Cancelation of gauge anomalies is a basic requirement for constructing grand unified models, 
and the usual assumption is that anomalies must cancel among spin | fermion fields. However, 
a 1985 paper of Marcus [1] showed that in principle an SU (8) gauge theory can be constructed 
with spin | Rarita-Schwinger fermions playing a role in anomaly cancelation, and we have recently 
constructed [2,|] a family unification model incorporating this observation. Using gauged spin | 
fields in a grand unification model raises the question of whether such fields admit a consistent 
quantum, or even classical theory. It is well known, from papers of Johnson and Sudarshan Qj 
and Velo and Zwanziger [4], and much subsequent literature (see e.g. Deser and Waldron [5J]), that 
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theories of massive gauged Rarita-Schwinger fields have serious problems. Does setting the fermion 
mass to zero eliminate these difficulties? 

The lesson we have learned from the success of the Standard Model is that fundamental fermion 
masses lead to problems and are to be avoided; all mass is generated by spontaneous symmetry 
breaking, either through coupling to the Higgs boson or through the formation of chiral symmetry 
breaking fermion condensates. So from a modern point of view, the Rarita-Schwinger theory 
with an explicit mass term is suspect. Several hints that the behavior of the massless theory 
may be satisfactory are already apparent from a study of the zero mass limit of formulas in the 
Velo-Zwanziger paper. First, in their demonstration of superluminal signaling, the problematic 
sign change that they find for large B fields (Eq. (2.15) of 4]) is not present when the mass 
is set to zero. Second, when the mass is zero, the secondary constraint that they derive (Eq. 
(2.10) of 4|]) appears as a factor in the change in the action under a Rarita-Schwinger gauging 
Sipfj, = D^e, with D^ the usual gauge covariant derivative. (This statement is not in [4], but is 
an easy calculation from their Eqs. (2.1)—(2.3), with the D^ of this paper their Hence, 

the constrained action in the massless gauged Rarita-Schwinger theory has a fermionic gauge 
invariance that is the natural generalization of the fermionic gauge invariance of the massless free 


Rarita-Schwinger theory. Third, their formula for the anticommutator (Eq. (4.12) of Bo in the 
zero mass case develops an apparent singularity in the limit of vanishing gauge field B, and so their 
quantization does not limit to the standard free theory quantization. This should not be surprising: 
since the massive theory does not have a fermionic gauge invariance, Ref. [4] does not include a 
gauge-fixing term analogous to that used in the massless case, but gauge fixing is needed to get 
a consistent quantum theory for a free massless Rarita-Schwinger held. So these observations, 
following from the equations in J 4 ], suggest that a study of the massless Rarita-Schwinger held 
coupled to spin-1 gauge helds is in order. 

In a different setting, massless Rarita-Schwinger helds appear consistently coupled to gravity 


as the gravitinos of supergravity, as discussed by Das and Freedman [6j]. Grisaru, Pendleton and 


van Nieuwenhuizen [7] have shown that soft spin | fermions must be coupled to gravity as in 
supergravity, in an analysis based on the free particle external line pole structure of spin § helds 
that do not have spin 1 gauge couplings. Their result does not rule out the possibility of gauged 
spin | helds with non-Abelian gauge couplings, since such fermion helds will in general not have free 
particle external lines. For example, a gravitino coupled in the fundamental representation of an 
unbroken non-Abelian gauge group will be conhned, and will not have the kind of external line pole 
structure used in the argument of [7]. Thus, the known connection of ungauged Rarita-Schwinger 
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fields to supergravity does not argue against the possibility that there could be a consistent theory of 
massless, gauged Rarita-Schwinger fields, so again a detailed study of this possibility is warranted. 

B. Outline of the paper, and summary 


With these motivations and background in mind, we embark in this paper on a detailed study 
of the classical and quantum theory of a minimally gauged massless Rarita-Schwinger field. In 
Sec. 2, we give the Lorentz covariant Lagrangian for a gauged four-component Rarita-Schwinger 
spinor field, derive the source current for the gauge field, and check that it is gauge-covariantly 
conserved. We also give the Lorentz covariant form of the constraints and of the fermionic gauge 
invariance, and of the symmetric stress-energy tensor. In Sec. 3 we argue that the massless, gauged 
Rarita-Schwinger theory exhibits a generalized fermionic gauge invariance, in the sense that the 
action principle does not uniquely determine the system time evolution. Since in the massless case 
left chiral and right chiral components of the field decouple, in Sec. 4 we rewrite the Lagrangian 
for left chiral components in terms of two-component spinors and Pauli matrices, which simplifies 
the subsequent analysis. We then give the Euler-Lagrange equations in two-component form, and 
use them to analyze the structure of constraints and the fermionic gauge invariance of the action. 
In Sec. 5 we specialize to the case of an Abelian gauge field (as in [4(]) and analyze the wavefront 
structure, showing that wave modes propagate with subluminal velocities; an extension of this 
discussion is given in Appendix B. In Sec. 6 we return to the general case of non-Abelian gauge 
fields. We introduce canonical momenta for the Rarita-Schwinger field components, which are used 
to define classical Poisson brackets, and discuss the role of the constraints as generators of gauge 
transformations under the bracket operation. We show that the constraints group into two sets of 
four, within each of which there are vanishing Poisson brackets. This permits application of the 
Faddeev-Popov method for path integral quantization, which we carry through in detail in Sec. 7 in 
\ko = 0 gauge. In Sec. 8 we give the Hamiltonian form of the equations of motion and constraints, 
and introduce the Dirac bracket. When a gauge fixing constraint is omitted, the Dirac bracket that 


we compute agrees with the anticommutator calculated in [3| and 4]. In Sec. 9 we study the Dirac 
bracket in its classical and quantum forms, and show in covariant radiation gauge that the quantum 
Dirac bracket has the requisite positivity properties to be an anticommutator; related details are 
given in Appendix C. In Sec. 10 we study the behavior under Lorentz transformations of covariant 
radiation gauge, of \ko = 0 gauge, and of the Dirac bracket and anticommutator. In Sec. 11 we 
analyze Rarita-Schwinger fermion scattering from a short range Abelian gauge potential, using an 
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analog of the distorted wave Born approimation. We show consistency of relativistic covariance 
with helicity state counting, but there is no small-coupling expansion giving an unmodified Born 
approximation for this process. In Sec. 12 we discuss areas for extensions of our results, and in 
Appendix A we summarize our notational conventions and some useful identities. We suggest that 
the reader skim through Appendix A before going on to Sec. 2, since things stated in Appendix A 
are not repeated in the body of the paper. 

Our conclusion from this analysis is that one can consistently gauge a massless Rarita-Schwinger 
field, at both the classical and quantum levels. This opens the possibility of using gauged Rarita- 
Schwinger fields as part of the anomaly cancelation mechanism in grand unified models, with 
anomalies of the spin ^ fields canceling against the spin | anomaly. 

II. LAGRANGIAN AND COVARIANT CURRENT CONSERVATION IN 

FOUR-COMPONENT FORM 

The action for the massless Rarita-Schwinger theory is 

S'CVV) = \ J d4x ^na U R^ au , 

R , QU =ie^( n r/s(D^) U , 

=d v ^ + , 

A u — A A t u . 

■^■vv u Av 5 

(i) 

with ^ au a four-vector four-component spinor, with four-vector index /j = 0,..., 3, spinor index 
a = 1, ...,4, and SU(n ) internal symmetry index u = 1, ...,n, with SU(n ) gauge generators t a, A = 
1, ...,n 2 — 1. Taking u to range from 1 to n means that, for definiteness, we are assuming that the 
spinors transform according to the fundamental representation of the SU(n ) internal symmetry 
group, but other representations and other compact Lie groups can be accommodated by assigning 
the internal indices u and A the appropriate range. Note that ia, A^ v , and D u all commute with 
the gamma matrices and the Pauli spin matrices from which the gamma matrices are constructed, 
and for an Abelian internal symmetry group, the indices u and A are not needed. Using 

VW = ^lpu i( CI 0 fa > ( 2 ) 

together with the adjoint convention (x\x 2 )^ = X 2 X 1 f° r Grassmann variables xii X 2 , it is easy to 
verify that S is self-adjoint. 
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From here on we will usually not indicate the spinor indices a, f5 and internal symmetry indices 
u, v explicitly, but they are implicit in all formulas. Varying S with respect to the Rarita-Schwinger 
fields, we get the equations of motion 

e^du^pj v =ge ppup ^ p A^t A lr 1 , 

= - ge^ p lv A^ p . 

(3) 


Re-expressed in terms of the covariant derivative, these are 

e^ p t„ lri = 0 , 

e ppvp lr] D v ^ p = 0 . 

(4) 


The g = 0 component of these equations gives the primary constraints 

e en > r K 7e =0 , 

e enr le D n ^ r =0 , 

(5) 

with e,n, r summed from 1 to 3. Contracting the equation of motion for if) with p and the 
equation of motion for ipp with g~ 1 D p , we get the secondary constraints 

e pp ^ p F pulv = 0 , 

e MU V^ p =0 , 

( 6 ) 


where we have introduced the gauge field strength 

F {W =g~ 1 [D l _ l ,D u \ = g~ l \F> p ,‘l5 u \ 
=dfj,A u d u A p + g[A p , A u ] , 


(7) 


which with the adjoint representation index A indicated explicitly reads 


F it = d u A u - d„A* + gfABcA^A^ 


(8) 
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Under a Rarita-Schwinger gauge transformation (with e a four-component spinor), which is a 
natural gauge field generalization of the fermionic gauge invariance for a free, massless Rarita- 


Schwinger field discussed in 


8 ], 

Vv 

Vv 


hI’h + SgiI’h , 


= D^e 
= gDu 


the action of Eq. © changes according to 


1 


SGStyJ = ~~ig / d 4 x 675 ( e mvp 'tnF l j V i}p ) + 


756 


(9) 


( 10 ) 


The factors bracketed in large parentheses are identical to the secondary constraints of Eq. ([ 6 ]) . 
We will argue in the next section that this implies that, even when coupled to gauge fields, the 
Rarita-Schwinger theory has a generalized form of fermionic gauge invariance. 

Adding the gauge field action 


s(K) = 4 


d A xF±F Atxv 


fJ,lS 


( 11 ) 


and varying the sum 5(4^) + S(A A ) with respect to the gauge potential, we get the gauge field 
equation of motion 

D v F a ^ = 8 V F A ^ + gf AB cA*F c ^ = gJ Afl , 


jAil =7'*/4^ i ' w 757^ j 4?/v 


( 12 ) 


A straightforward calculation using Eqs. ([3|) shows that the gauge field source current J A ^ obeys 
the covariant conservation equation 

+ gf ABC Ap c » = 0 , (13) 

as required for consistency of Eq. (fl2]l . So from the Rarita-Schwinger and gauge field actions, we 
have obtained a formally consistent set of equations of motion. 

In addition to the gauge field source current, there is an additional current that obeys an 
ordinary conservation equation, 




757r?^ 


P ’ 


=0 


( 14 ) 
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In the massive spinor case, Velo and Zwanziger [4] argue that the analogous current, within the 
constraint subspace of Eq. (J5|) , should have a positive time component. In the massless case 
we see no reason for this requirement, since Eq. (114(1 is the fermion number current, and its time 
component, giving the fermion number density, can have either sign. However, we shall use parts of 


the positivity argument of 


4| later on in discussing positivity of the Dirac bracket anticommutator. 


The symmetric stress-energy tensor for the free massless Rarita-Schwinger has been computed 


by Das 


(see also Allcock and Hall [10(]). Changing ordinary derivatives to gauge covariant 


derivatives, Das’s formula becomes 


* A [ivp 


T ar —_ c 

J RS - 4 e 

* ^ 


^a 75(7 t ^ + l a ^) D^ P 


(15) 


(This formula can be made manifestly self-adjoint by replacing D v by \{D V — *D U ), but this is not 
needed to verify stress-energy tensor conservation.) Adding the gauge field stress-energy tensor, 


rpCTT 

1 gauge 


= ~V ar F^F Ax ^ + f£°F AXt , 


(16) 


a lengthy calculation, using Eq. (1131) together with identities and alternative forms of the equations 
of motion given in Appendix A, shows that the total tensor is conserved, 


<9(7 (?RS + ^gause) — 0 


(17) 


III. GENERALIZED GAUGE INVARIANCE OF THE RARITA-SCHWINGER ACTION 


In the most familiar gauge invariant theories, such as Abelian or non-Abelian gauge fields, the 
Lagrangian density is invariant under a gauge transformation on the fields. These theories exhibit 
what one could term “strong” gauge invariance. In a weaker form of gauge invariance, which occurs 
for the free Rarita-Schwinger equation, the Lagrangian density changes by a total derivative under 
a gauge transformation of the fields, and so only the action is gauge invariant. Characteristic 
features of this case have been studied by Das [9j|. We argue in this section that there is a still 
weaker form of gauge invariance, obeyed by the massless Rarita-Schwinger equation with Abelian 
or non-Abelian gauging, in which under a gauge transformation the Lagrangian changes by a total 
derivative plus terms which vanish when initial value constraints are obeyed. 

In his seminal analysis of constrained systems, Dirac [11] classifies as “first class” constraints 
the maximal set of constraints that have vanishing mutual Poisson brackets, and notes that “Each 
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of them thus leads to an arbitrary function of the time in the general solution of the equations of 
motion with given initial conditions”. Elaborating on this, he notes that “Different solutions of the 
equations of motion, obtained by different choices of the arbitrary functions of the time with given 
initial conditions, should be looked upon as all corresponding to the same physical state of motion, 
described in various way (sic) by different choices of some mathematical variables that are not of 
physical significance (e.g. by different choices of the gauge in electrodynamics or of the co-ordinate 
system in a relativistic theory.)” 

These remarks suggest that gauge invariance, in its most general form, corresponds to an ar¬ 
bitrariness in the time evolution of a system, in the sense that the future evolution of the system 
is not uniquely determined by the initial conditions and the Euler-Lagrange equations following 
from the action principle. Under this generalized definition, the Rarita-Schwinger equation with 
coupling to gauge fields has a fermionic gauge invariance. To see this, we note the Euler-Lagrange 
equations yield equations of two types. The first are the time evolution equations contained in Eq. 
([3]), that determine the field variables at time t + At from those initially given at time i. Since we 
will see that ip$ is related by a constraint to ip, the evolution equations take the form 


ip(t + At) = ip(t) + At^ ip(t),ip 0 (t) + 0((At) 2 ) 


(18) 


The second are the primary and secondary constraints of Eqs. ([5]) and ([6]), which constrain the 
initial field values at time t. If we make the gauge transformation of Eq. Q at time t, with 
gauge parameter e of order At, we have seen that the action at time t changes, to first order in e, 
according to Eq. m- So when the constraints at time t are applied, the change in the action is 
0((At) 2 ^ . The effect of the gauge shift Script) of order At is to change the time evolution equation 
Eq. (fTHl) to 


ip(t + At) =ip(t) + S G ip(t) + At^ ip(t) + d G ip(t),ipo{t) + 5 G ipo(t) +0((At) 2 ) 
=ip(t) + S G ip(t) + At^ ip(t),ip 0 (t) +0((At) 2 ) . 


(19) 


In other words, the error induced in the Atdip/dt term by the gauge transformation is 0((At) 2 ). 
Hence the gauge transformed variables are equally good candidates for time evolution to t + At, 
leading to a state of the system at t + At that is no longer unique; for each independent gauge 
change at time t we get a different set of evolved variables at time t + At. One can now evolve 
the variables at t + At to time t + 2At, subject to the constraints at time t + At, allowing one to 
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make a further independent gauge transformation at time t + At, and thus introducing a further 
non-uniqueness into the time-evolved solution. Continuing in this fashion, forming a Riemann sum 
in the limit At —>• 0, the 0((At) 2 ) corrections drop out, and so there is no change in the action 
associated with the gauge-modified system evolution. Thus one can introduce arbitrary functions of 
time into the evolved solution, corresponding to each independent fermionic gauge transformation 
of Eq. ([9]). In order to get a unique time evolution path from the initial data at time t, one must 
impose a gauge fixing condition, that selects at each time one member out of the equivalence class 
of equal action configurations. 

Therefore, even though only the constrained action in the gauged Rarita-Schwinger theory is 
invariant under fermionic gauge transformations, we shall refer to this as gauge-invariance in the 
sense just described, and will impose a gauge fixing condition to break this invariance. The gauge 
fixing condition will also serve to give the correct helicity counting for the Rarita-Schwinger fields, as 
we shall see in detail below. We shall also discuss below how the first class constraints, which follow 
from repeatedly differentiating the original primary constraints with respect to time, play a role 
as gauge transformation generators, in agreement with what is found in strongly gauge invariant 
systems. However, we shall see that a novel feature emerges in the Rarita-Schwinger theory: When 
one insists that for each first class fermionic constraint there should be a corresponding adjoint 
constraint, doubling the constraint count of the theory, then the first class constraints become 
second class constraints in the Dirac classification. 

IV. LAGRANGIAN ANALYSIS FOR LEFT CHIRAL SPINORS IN TWO-COMPONENT 

FORM 

Although we could continue with the four-component formalism to study constraints, the Hamil¬ 
tonian formalism, and quantization, it will be more convenient to first reduce the four component 
equation to decoupled equations for left and right chiral components of ^ (with a the spinor 
index and with the internal symmetry index implicit). Since these are related by symmetry, we 
can then focus our analysis on the two-component equations for the left chiral component, which 
is the component conventionally used in formulating grand unified models (see, e.g. jsj]). 

We convert the action of Eq. © to two-component form for the left chiral components of 
using the Dirac matrices given in Eqs. (jA2l) and (1A41) . Defining the two-component four vector 
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spinor and its adjoint by 

PiM’Z = ( ^ ) , M = 0,1,2,3, a = 1,2 , 

^U PL= iyU o) , 


( 20 ) 


the action decomposes into uncoupled left and right chiral parts. The left chiral part, with spinor 
indices a suppressed, is given by 

S(^) = i J d 4 x[-Tj<T • D x T + # • a x D^ 0 + T t -ZlxT-T t -<7X D 0 $] . (21) 

Varying with respect to we get the Euler-Lagrange equation 

0 = V = a x o + D x T - a x D 0 $ , (22) 

while varying with respect to Tg we get the primary constraint (given in four-component form in 
Eq. ©) 

D = X = a • z) x T . (23) 

A second primary constraint follows from the fact that the action has no dependence on d^/dt, 
which implies that the momentum conjugate to Tq vanishes identically, 

P* t = 0 • (24) 

Contracting V with a and with g _1 D, and using the covariant derivative relations of Eq. (1A14I) . 
we get respectively 


a ■ V =2 id + X , 


g 1 D -V =iu + g 1 D 0 x , 


(25) 


with 


9 =a ■ .DTo — Dqo ■ T , 
oj =a ■ o — (B + a x E) • T 


(26) 
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Since the Euler-Lagrange equations imply that V and x vanish for all times, we learn that 6 and uo 
vanish also for all times. Since 9 involves a time derivative, its vanishing is just one component of 
the equation of motion for . But io involves no time derivatives, so it is a secondary constraint 
that relates d'o to (given in four-component form in Eq. ([6])). For each of the above equations, 
there is a corresponding relation for the adjoint quantity. 

The equation of motion V = 0 can be written in a simpler form by using the identities of Eqs. 
(jAlOl) and (IA11I) as follows. Using Eq. ()A10f) to simplify 0 = <7 x V — iV, we get an equation for 

A)*, 

D 0 $ = D\ h 0 + ~ [— a x (D x T) + iD x \k] . (27) 

A further simplification can be achieved by incorporating the primary constraint x = 0, through 
applying Eq. (1A11D to A = D x T, 

0 = ax = vv-(.Dx$) = Dx$-iax(Dx$) . (28) 

Using this to replace the first term in square brackets in Eq. (1271) we get the alternative form of 
the equation of motion, valid when the constraint x = 0 is satisfied, 

£) 0 T = .DTo + iD x T . (29) 

Writing the gauge field interaction terms in Eq. (1211) in the form 

Sint(’JV) = f / d A x(A*J m + A b • J B ) , (30) 

we find the left chiral contribution to the current of Eq. (1121) in the form 

jAO = _ . g x ^ s 

J A =^ 0 tAd x T + x at a 'ho - x tA$ ■ 

(31) 


Replacing t,A by —i, we find the corresponding singlet current in the form 

J° =*Tt • a x T , 

J = - i(Tj(T xl + l*x <tT 0 - x f ) . 

(32) 


For the energy integral computed from the left chiral part of the the stress-energy tensor of Eq. 
(fl5l) . we find 


d 3 xT™ = - 


1 


H = — 


d 3 x& ■ D x T 


(33) 
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To conclude this section, we verify that the action of Eq. (1211) has a fermionic gauge invariance 
on the constraint surface oj = 0 , = 0, as already seen in covariant form following Eq. ([9]). 

Letting e be a general space and time dependent two-component spinor, we introduce the fermionic 
gauge changes 

T —+ <5 G T , <5 G T = De , 

*o —+ ‘Sg'I'o , feTo = D 0 e , 

(34) 

and their adjoints, which are the left chiral form of the gauge change of Eq. @. Substituting this 
into Eq. (ED), integrating by parts where needed, and using Eqs. (IA14I) to simplify commutators 
of covariant derivatives, we find that Eq. (HOD takes the two-component spinor form 

dcSi'S'fj,) = ^ig j d A x{u*e- e^u) . (35) 

Hence the action on the constraint surface uj = = 0 has a fermionic gauge invariance. Another 

gauge invariant, on the constraint surface y = y^ = 0, is the fermion number, given by the space 
integral of the time component of the singlet current of Eq. (|32D . f d 3 xJ° , which has the gauge 
variation 



However, neither the equation of motion, the constraints y and w, the non-Abelian “charge” 
f d 3 xJ B0 , nor the integrated Hamiltonian H are gauge invariant in the interacting case. Using <5 G 
to denote gauge variations, we have 

5qV = — ig(B + if x E)e , 

5q0 = — igcr ■ Ee , 
fcx = -igo - Be , 

Sqoj =cr • BD$e — {B + a X E) ■ De , 

Sg j d 3 xJ B0 =g j d 3 x(e^[A, ts] ■ a x T + 4^ x a ■ [tB,A]e^j , 

5 g H =\ig J d 3 x{^ ] • Be - e ] B ■ T) . 

(37) 

The last line shows that in the interacting case, there is no total energy integral that is invariant 
under the fermionic gauge transformation, much as there is no generally covariant gravitational 
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total energy. The only global fermionic gauge invariants are the action integral, and the fermion 
number integral. 

To break the gauge invariance we must introduce an additional constraint, in the form 

/W = 0 , (38) 

with / a scalar function of its argument, such as f = D ■ T (a gauge covariant radiation gauge 
analog) or / = (B + a x E) ■ T (which by Eq. (1261) corresponds, when a ■ B is invertible, to 
To = 0.). This constraint, together with the x constraint, leaves one independent two-component 
spinor of the original three in T, corresponding to the physical massless Rarita-Schwinger modes 
propagating in the gauge held background. We will limit ourselves to considering linear constraints 
of the general form 

/ = L ■ T , (39) 

so that for gauge covariant radiation gauge we have L = D, and for To = 0 gauge we have 
L = B + a x E. Both of these choices of L play a special role in our analysis, so we shall examine 
them in more detail. 

We consider first the gauge covariant radiation gauge. We note that since 

a ■ Da • T = D ■ T + i\, (40) 

the primary constraint X = 0 implies that 

a ■ Da • T = D ■ T . (41) 

Hence when a ■ D is invertible, which is expected in a perturbation expansion in the gauge coupling 
g, the covariant radiation gauge constraint D ■ T = 0 implies that 

a ■ T = 0 . (42) 

Conversely, Eqs. m and (1411) show that D ■ T = 0 and a ■ T = 0 together imply the primary 
constraint x = 0, and also a ■ T = 0 and x = 0 imply D ■ T = 0. 

We next note that on a given initial time slice, covariant radiation gauge is attainable. Under 
the gauge transformation of Eq. (1M|) , we see that 


D ■ T -► D ■ T + {Dfe 


(43) 
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Hence when ( D ) 2 is invertible, which we expect to be true in a perturbative sense, then we can 
invert (D) 2 e = —D • 4', to find a gauge function e that brings a general T to covariant radiation 
gauge. Since 

(a • Df = (D) 2 + go ■ B , (44) 

the conditions for a ■ D to be invertible, and for ( D ) 2 to be invertible, are related. For generic 
non-Abelian gauge fields both of these operators should be invertible, but there will be isolated 
gauge field configurations for which a ■ D has zeros. 

However, although covariant radiation gauge can be imposed on any time slice, it is not preserved 
by the equation of motion for 4'. To see this, let us consider the simplified case in which the gauge 
potential is specialized to Ho = 0 and do A = 0, so that only a static B held is present. From Eq. 
(l29l) we have 

d 0 (D ■$) = (D) 2 'I> 0 +gB ■$ = ((D) 2 (a ■ B)- 1 + g)B-^ /0 . (45) 


Hence at each infinitesimal time step, we must make a further infinitesimal gauge transformation 
to maintain the gauge covariant radiation gauge condition. This means that we cannot use the 
9 = 0 time development equation 


<j • ZJTo = Dqij • T 


(46) 


together with Eq. (142(1 . and the assumption that a ■ D is invertible, to conclude that covariant 
radiation gauge also implies that T 0 = 0. 

We can, however, impose To = 0 as an alternative gauge condition. Under the gauge transfor¬ 
mation of Eq. (j34j) , To can be reduced to zero by taking 

( 


e = — 


V(t,x) 1 f duV(u, x)To(tt, x) , 
Jo 


U(f,x)= expig / duAo(u,x) 


(47) 


Once To has been gauged to zero, we can conclude from Eq. (|26l) that (B + a x E) ■ T = 0. 
Just as we found for covariant radiation gauge, To = 0 gauge is not preserved by the equation of 
motion for T. To see this, let us again consider the simplified case in which the gauge potential is 
specialized to Hq = 0 and do A = 0, so that only a static B held is present. From Eq. (1291) we have 


d 0 {B- T) =iB -D x T /0 


(48) 
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So again, at each infinitesimal time step, we must make a further infinitesimal gauge transformation 
to maintain the gauge condition 'ho = 0. 

To summarize this discussion of gauge fixing, we see that when the gauge fields are nonzero, we 
can have either ITT = 0 or To = 0, but not both. In the free field case, these two conditions can be 
imposed simultaneously and are preserved by the time evolution of T, so we recover the constraints 


V • ip = 0, a ■ ip = 0, and i/jq = 0 used in the discussion of 


8| for the free Rarita-Schwinger case, 


V. PROPAGATION OF A RARITA-SCHWINGER FIELD IN AN EXTERNAL 
ABELIAN GAUGE FIELD: ABSENCE OF SUPERLUMINAL PROPAGATION 


We specialize now to the case of a Rarita-Schwinger spinor propagating in an external Abelian 
gauge field, as studied by Velo and Zwanziger [4], For an Abelian gauge field, 

1 a ■ B 


(49) 

a-B (B ) 2 

and so o ■ B is invertible as long as ( B ) 2 > 0, which we assume. Provided the Lorentz invariant 
expression (B) 2 — (E) 2 is positive, ( B) 2 will be positive in any Lorentz frame. In discussing 
undamped wave propagation we will not use the inequality ( B) 2 — ( E ) > 0, but in treating damped 
longitudinal mode propagation in Appendix B, we will assume that ( E) 2 /(B) 2 is small, as motivated 
by the fact that when ( E) 2 is of order ( B) 2 the vacuum is highly unstable against pair creation. 
(Strictly speaking, the vacuum is stable against pair production only when E ■ B = 0 and ( B) 2 — 
( E ) 2 > 0, that is, when there is a Lorentz frame in which the Abelian field has vanishing E [12j].) 
Given that ( B ) 2 > 0, we can solve the constraint u = 0 of Eq. (1261) for To, giving 

Q • T 


T n = 


{BY 


(50) 


where we have defined 


Q = a • B(B + axE) = BxE + Ba ■ (B + iE ) — iB ■ Ea . (51) 

Substituting the solution for To into Eq. P§1) . we get an equation of motion for T by itself, 

(52) 


- -*Q ■ T 

D 0 T = D __ IiDxT 


(BY 

To determine the wave propagation velocity in the neighborhood of a point a?*, we need to 
calculate the equation for the wavefronts, or characteristics, at that point. Writing the first order 
Eq. (1521) in the form 

■ r Q* ■ $ 


<9 0 T = V- 


(B* 


+ i\7 x T + A[T,x*,x] , 


(53) 
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with B * and Q* the values of the respective quantities at x*, we see that A[\Il,x*,xl involves no 

. n n 

first derivatives of T at f*, and so is not needed [13], [14|] for determining the wavefronts of Eq. 
(1291) . The reason is that when taking an infinitesimal line integral of Eq. (1531) , according to 

lim / [d 0 T = ...] , (54) 

J-5 

discontinuities across wavefronts contribute through the first derivative terms, but when the exter¬ 
nal fields are smooth the term A[T,x*,x] makes a vanishing contribution as 5 —> 0. Dropping A, 
and multiplying through by ( B *) 2 , we get the equation determining the wavefronts in the form 

(S*) 2 d 0 T = VQ* • T + i{B*fV x T . (55) 

By similar reasoning, the constraint % can be simplified, for purposes of determining the wavefronts, 
by replacing D by V, giving 

0 = e? • V xf , (56) 

and the covariant radiation gauge condition similarly simplifies to 

0 = V • T . (57) 

Since these are now linear equations with constant coefficients, the solutions are plane waves, 
and without loss of generality we can take the negative z = x 3 axis as the direction of wave 
propagation. So making the Ansatz 

T = C ex.p(iQt + iKz) , (58) 

Eq. (1551) for the wavefronts or characteristics takes the form 

0 = F= (B*) 2 nC - KzQ * • C - i(B*fKz x C , (59) 

with z a unit vector along the z axis. Similarly, the constraint Eq. (1561) becomes an admissability 
condition on (7, 

0 = a ■ z x C , (60) 

and the gauge fixing condition V ■ T = 0 imposes the further condition on C 

0 = z ■ C , (61) 

but for the moment we will analyze the solutions of Eq. (159[) without assuming this additional 


constraint. 
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Writing F m as a matrix times C n (and dropping the subscripts *, which are implicit from here 
on) we have 


Fm, — -\ r rn n Cn 


Nmn — (B) B5 m 3Qn 


i(B') Ke m 3 n 


( 62 ) 


The equation for the characteristics is now 


det(lV) = 0 , (63) 

since this is the condition for Eq. (|59l) to have a solution with nonzero C. However, since evalua¬ 
tion of the determinant shows that it factorizes into blocks that determine C\ 2 and a block that 
determines C 3 , a simpler way to proceed is to work directly from the equations F m = 0 , which 
decouple in a corresponding way. Calculating from Eq. 11591) . we find 

0 =Fl’ 1 = (B) 2 (HC'J'T + iKC\ A ) , 

0 = F 2 t ’ 4 ' = (B) 2 (yLCl A - iKC\ A ) , 

0 =f} 1 = (B) 2 nc} 1 -K(Q-C)^ 1 , 

(64) 

where "f, 4- indicate the up and down spinor components, labeled in Eq. (I2Q[) by a = 1, 2. Similarly, 
the constraint Eq. (f60l) becomes 0 = — CT 1 C 2 + ogCi, that is 

C$=iCl , 

ni — __ Ani 
L / 2 — , 

(65) 

-f- I 

with no corresponding condition on C 3 . 

The first two lines of Eq. (1641) together with Eq. (165(1 have the solution 

C\ =C , cl = iC , n = K , 
cl =c, c'l = -iC , n = —K , 

( 66 ) 


with C arbitrary, corresponding to waves with velocity of magnitude 112/77) = 1 . The effect of 
a gauge change T —> ^ + De, e = (E/(zA')) exp(zllf + iKz ) on the characteristics is to shift 
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Cg"*" —> + thus replacing the third line of Eq. (|64l) by the two-component spinor equation 

((i?) 2 0— KQs) (Cs+E)—K(Q\C\+Q 2 C 2 ) = 0. If we choose E to cancel the term A" (QiCi+l^C^)) 
that is, 

e = , (67) 

{B) 2 n - kq 3 

then this equation is solved by C 3 = 0 (corresponding to enforcing the covariant radiation gauge 
conditions V • T = <7 ■ \E' = 0), and the demonstration that there is no superluminal propagation is 
complete. In Appendix B, we continue this discussion without assuming a specific gauge condition, 
and show that there are still no propagating modes with superluminal velocities. That is, the 
purely longitudinal mode with C\ = C 2 = 0, C 3 ^ 0, which is eliminated by the gauge fixing 
condition, does not propagate superluminally. 


VI. CANONICAL MOMENTA, CLASSICAL BRACKETS, AND GAUGE 

GENERATORS 


We return now to the general non-Abelian gauge field case, and introduce the canonical mo¬ 
mentum conjugate to T, defined by 


5 dLS 1 ft - 
P = -— = -T T x a 

9(9 0 T) 2 

which can be solved for Tt using the final line of Eq. (IA11D . 


( 68 ) 


#t = iP-P 


X tT 


(69) 


We will use Eq. (|69D when computing classical brackets involving 4 /t using the formula of Eq. 
(1A17D . Eq. (1681) can be written as an explicit matrix relation for the six components of P and T', 


/pa 


( 0 

0 

1 

0 

0 

-z\ 


( 

pi 
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0 

0 

-1 
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0 
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0 

-1 
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V -* 

0 

-1 

0 

0 

0 ) 




(70) 


showing that they are related by an anti-self-adjoint matrix with determinant —1/16. 
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The four constraints introduced in Sec. 4 are 


=P fa ’ 

<f>2 =(<? • Bfa 1 oj = — (<t • B)~ l (B + a x E) ■ ^ , 

fa =X = & ■ D x , 

(j>4 =L ■ T 

(71) 

In writing these we are assuming that a • B is invertible in the non-Abelian case. We are writing 
the gauge fixing condition as a general linear gauge fixing constraint L ■ T, so as to keep track of 
which terms in the final answers arise from gauge fixing, which is not evident if we specialize by 
replacing L by D or B + a x E. The constraints of Eq. (1711) . including the gauge fixing constraint 
fa, are all first class in the Dirac classification, since they have vanishing mutual classical brackets. 
This is a consequence of the fact that starting with a constraint depending on T but not on 4/1, 
and taking an arbitrary number of time derivatives, one still has a constraint depending only on 
T. 

To preserve the adjoint properties of the Rarita-Schwinger equation, for each of these four 
constraints we must impose a corresponding adjoint constraint. Using Eq. (|69l) to express in 
terms of P, we write these as 


Xi =C p *t) t = -P* o , 

X 2 =w t (a-R)“ 1 = Tj -P- [i(B + a x E) - & x (B + a x E)](a ■ B)~ l , 

X3 =X f = 2 P-t) , 

X 4 ■ L = P ■ {iL — a x L) 

(72) 


(The reason for the minus sign in the definition P ^t = — p|, q will be given in Sec. 8 where we 
discuss the Hamiltonian form of the equations.) The constraints 4> a are implicitly 2 n component 
column vectors, and the adjoint constraints Xa are implicitly 2 n component row vectors, with 2 n 
arising from the product of a factor of 2 for the two implicit spinor indices, and a factor of n for 
the n implicit SU(n) internal symmetry indices. 

When L = D, we see that <p4 becomes $4 = D ■ T, and x ,4 becomes %4 = • ^4 — P • cj x *5 = 

(t/2)x3 -E-ffxB. So a special feature of covariant radiation gauge, which will be exploited later, 
is that the constraints fa, fa are contractions ofiJxD and D with T, and the constraints X 3 , X4 
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are contractions of linear combinations of the duals D and a x D with P. That is, in covariant 
radiation gauge the constraint spaces selected by X3>X4 and are similar. 

We can now compute the classical brackets of the constraints. We see that the brackets of the 
4>s and x s vanish among themselves, 


\4>aAb\c =0 * 

[Xa,Xb\c = 0 , 
a, b =1,..., 4 


(73) 


On the other hand, the brackets of the cj)s with the x s give a nontrivial matrix of brackets M, 
which has a nonvanishing determinant, 


M a b{x, y) =[4> a (x),Xb{y)]c + 0 , 
det M . 


(74) 


Thus, in terms of the Dirac classification, the original first class constraints (j) a have become second 
class, not from adding new constraints that follow from differentiation with respect to time or from 
imposing gauge fixing conditions, but rather from adjoining the adjoint set of constraints. This is 
a feature of the Rarita-Schwinger constrained fermion system that has no analog in the familiar 
constrained boson systems such as gauge fields. 

Evaluating the brackets shows that M has the general form 

/0 -1 0 0 ^ 

1 U S T 

M = 

0 V A B 
\ 0 W C V j 

where in the SU(n) gauge field case, each entry in M is a 2 n x 2 n matrix (corresponding to the 
fact that 4> a is implicitly a 2 n component column vector, and Xb is implicitly a 2 n component row 
vector). Evaluating det M by a cofactor expansion with respect to the elements of the two unit 
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matrices ±1, we see that the submatrices U , S , T, V, W do not contribute, and we have 


detM =detJV 



(76) 


So we need to only evaluate the brackets M 33 = A, M 34 = B , M 43 = C, M 44 = P, giving 

A = — 2iga ■ B(x)5 3 (x — y) , 

B = -2 D g - L s 6 3 (x-y) , 

C =2L$ ■ D^5 3 (x — y) , 

V = (z(4) 2 + 3 • (4 x 4))(s - V) ■ 

(77) 


When L = D, these become 


A = — 2iga ■ B(x)5 3 (x — y) , 

B = -2(D s )H 3 {x-y) , 

C=2(4) 2 5 3 (x-y) , 

V =i((4) 2 - 5 "? • B(x))5 3 (x - y) 


(78) 


Reflecting the fact that the (f> a and Xa are adjoints of one another, together with the fact that 
the matrix relating \IR to P is anti-self-adjoint (see Eq. dZQD), these matrix elements obey the 
adjoint relations 

M a b(x,y)^ = -M ba (y,x) . (79) 


Applications of these bracket and determinant calculations will be made in the next two sections. 

To conclude this section, we note that the constraints y, 4 Pq, 0 , P xj/ t play the role of gauge 
transformation generators. For example, we have (with common time argument t suppressed) 


d 3 x-x\x)e(x),^(y) 


=Dye(y) , 


J C 


- J d 3 xPq, o (x)D 0 e(x),^ 0 (y) 


c 


=D 0 ffe(y) 


(80) 


So the fermionic gauge transformation is a canonical transformation. 
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VII. PATH INTEGRAL QUANTIZATION IN 'Lq = 0 GAUGE 


In studying quantization, we will specialize to the case where the external gauge potentials, 
and hence D, are time independent, since the simplest discussions of constrained systems assume 
time-independent constraints. This assumption can be dropped when the gauge field is quantized 
along with the Rarita-Schwinger field, but a more complex system of constraints and constraint 
brackets will then appear; we defer this extension to a future investigation. 

When the constraints are time independent, the classical brackets of Eqs. (1731) and (1741) have 
the form needed to apply the Faddeev-Popov [151 ] method for path integral quantization. (This has 
been applied in the free Rarita-Schwinger case by Das and Freedman [8] and by Senjanovic j3].) 


The general formula of [15l] for the in to out S matrix element (up to a constant proportionality 
factor) reads 


(out|S|in) oc f exp ( iS(q,p )) n 
d t 

dp{q,p) = n 5(Xa)5(</>a)(det[(/) a ,Xfe]) 5 n dpidqi 


(81) 


where £ = 1 when all canonical variables are bosonic, and £ = — 1 in our case in which all canonical 
variables are fermionic, or Grassmann odd. In applying this formula, we note that since the action 
S of Eq. (12T1) and the bracket matrix M of Eqs. ([7i]) - ([77]) are independent of P^ 0 and P ^\, we can 
immediately integrate out the delta functions in these two constraints. Also, since the canonical 
momentum P is related to \IR by the constant numerical transformation of Eqs. (1691) and (1701) . we 
can take IIP as the integration variable instead of P. up to an overall proportionality constant. So 
we have the formula 

(out | S'| in) oc j exp (i— j d 4 x[— Tg(j • D x T + • a x DT o + DoT]) 


t,X 


with 


d/z(tf 0 , + ) = ) (det[<l> a ,Xb])- 1 d*od*ld$d& 


(82) 


(83) 


\c=2 


with dTo and dTj each a product over the spinor components, and dT and dT^ each a product 
over the spinor-vector components. 
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As our next step, we can carry out the integrations over To and Tq, using the delta functions 
5((f> 2 ) and <5(x2) ■ This leaves the formula 

(out|5|in) oc J exp (iX j d 4 x[—^ ■ (B + a x E){a ■ B)~ l & ■ D x T 

+# • (j x D(a • + 1? x £) • f + f f • -D x f - ^ • ff x A)'?]) 


t,X 


(84) 


with 

d/r(T,T t ) = ^Jj5(xc)<5(0 c )j (det[<)) a ,x fe ]) _ 1 dTdT t , (85) 

so that the only remaining constraints ^> 3 , 4 , X 3,4 are the ones used in constructing the determinant 
det [<t> a ,Xb]- 

We now confront a dilemma: If we impose the constraint x = 0 coming from d(x), the first 
term in the exponent drops out, and we are left with 

(out|5|in) oc J exp (i— J d A x^ • [a x D(a ■ + <7 xE)-T + Z)xT — <rx DoT]) 

xH^OM f ) , 

t,x 

( 86 ) 


in which the coefficient of \IT in the exponent is —V, and so at the stationary phase point of the 
exponent with respect to variations in \fd we get the correct equation of motion. On the other hand, 
if we integrate the first term in the exponent of Eq. (l 86 l) by parts, it becomes —(i/2) f d A x\\d ■ 
B) l {B + g x E) ■ T, which is set to zero by the constraint x^ = 0 coming from <5(x^)! Are we 
not allowed to integrate by parts here because the factor (a ■ B)~ A {B + a x E) is not continuously 
differentiable? This dilemma is avoided by working in Tq = 0 gauge, since as we have seen, then 
{B + a x E) ■ T = 0 and the troublesome term is absent from both the path integral and the 
equation of motion. In To = 0 gauge, we end up with the elegant formula 


(out | S | in) oc J exp (i- J d 4 xT^ ■ [D x T — a x I?o^]) 


x[]^(T,Tt) 

t,X 


(87) 
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In using this formula, the customary procedure [17f] would be to put the bracket matrix that is the 
argument of the determinant back into the exponent by introducing bosonic ghost fields cf)Q. But 
since for To = 0 gauge L = B + a x E, (jr G N(j)Q is only linear in derivatives of (j) G , so one would 
not get ghost propagators of the standard bosonic form. 


VIII. HAMILTONIAN FORM OF THE EQUATIONS AND THE DIRAC BRACKET 


An alternative way to quantize is to transform the Lagrangian equations to Hamiltonian form, 
and to take the constraints into account by replacing the classical brackets by Dirac brackets. In 
carrying this out, we will simplify the formulas by making the gauge choice Ao = 0 for the non- 
Abelian gauge fields. This gauge choice is always attainable, and leaves a residual non-Abelian 
gauge invariance A —> A + DA(x),with the gauge parameter A time independent. The Hamiltonian 
will then be covariant with respect to this restricted gauge transformation. For the moment, in 
discussing the canonical Hamiltonian and bracket formalism, we will allow A to be time dependent, 
so that E ^ 0. But when we turn to the Dirac bracket construction corresponding to a constrained 
Hamiltonian, we will assume a time-independent A as in the path integral discussion, corresponding 
in A 0 = 0 gauge to E = 0. (If we carry along the Ao term in the formulas, which we have done 
as a check, then time-independent fields would not require E = 0. So this specialization can be 
avoided at the price of somewhat lengthier equations.) 

From the action S(\ F^) = f dtL ( x of Eq. 1211) and the canonical momentum of Eq. (1681) . we 
find the canonical Hamiltonian to be 


H = 


d 3 xd 0 T P-L 


= ~ J d 3 x[- x &la ■ D x T + 'F t • a x DT 0 + • D x T] 

= — ^ J d 3 x[— Tju • D x T + ( iP — P x a) ■ (a x DTq + D x T)] 


( 88 ) 


where in the final line we have used Eq. (1691) to express in terms of P. 

We can now compute the classical brackets of various quantities with H. From 


dt 


= [T, H]c = ^[*(ct x DT o + DxT)-<7x(<7X DT 0 + D x T)] 
=DT 0 + ^[-a x (D x T) + iD x T] , 


(89) 
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we obtain the T equation of motion in the form given in Eq. (1271) . Similarly, from the bracket of 
P with H we find the equation of motion for Tb Turning to brackets of the constraints with H, 
starting with P T a , we find 


dP 9 

dt 


t 

0 





(90) 


and so P xJ a = 0 for all times implies that x = 0. For the total time derivative of X- we have 


^ + [x, H ]c = 3 x ■ T + [x,H] c = -igu , 


(91) 


and so x = 0 for all times implies that u defined in Eq. (126 jl vanishes. Since uo contains a term 
proportional to To, to continue this process by calculating the time derivative of u, we must obtain 
dTg/dt from a bracket of To with H (and similarly for dTj/di). This requires adding to H a term 


A H = - 


d 3 x 


P* o 


dTo 

dt 


+ P, 


K 


dt 


(92) 


Requiring AH to be self-adjoint then imposes the requirement 


4 0 p^ 0 


(93) 


which was noted following Eq. (1721) . We remark again that the chain of successive brackets with 
H starting from P r a and continuing to x,^,— leads only to constraints involving T and To but 
never their adjoints. The doubling of the set of constraints, which turns the first class constraints 
into second class ones, comes from requiring that the adjoint of each fermionic constraint also be 
a constraint, not from taking successive brackets with H. 

We are now ready to implement the Dirac bracket procedure. The basic idea is to change the 
canonical bracket [F, G]c to a modified bracket [F. G\d, which projects F and G onto the subspace 
obeying the constraints, so that the constraints are built into the brackets, or after quantization, 
into the canonical anticommutators. The constraints can then be “strongly” implemented in the 
Hamiltonian by setting terms proportional to the constraints to zero. In doing this we will choose 
the gauge To = Tj = 0, as we did in the path integral discussion. Since after integration by parts 
the second line of Eq. (1881) takes the form 

H = ~I d 3 x[-Tjx - X^o + ^ -DxT] , (94) 

this eliminates problems associated with the fact that To, Tq have as coefficients the constraints 
X f , X respectively. Setting constraint terms to zero in Eq. (1941) we see that the constrained 
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Hamiltonian is just 


H = - 



d 3 x& ■ D x $ 
d 3 x{iP — P x a) ■ D x T 


(95) 


which coincides with the energy integral computed in Eq. (ITill from the stress-energy tensor. 

We proceed now to calculate the Dirac bracket for the case when F = F(^) and G = G(T, \Ib); 
the case when F = F(^) can then be obtained by taking the adjoint, and the case when F = 
can be obtained by combining the extra bracket terms from both calculations. When 
F has no dependence on \Ib, it has vanishing brackets with the constraints cj) a of Eq. (1711) and 
nonvanishing brackets with the constraints Xa of Eq. (1721) . The Dirac bracket then has the form 
(see Eqs. (|A20I) and (lA21j) for why M -1 appears) 


[F,G\d = [F, G]c Xa]c M p b , G\ , (96) 

a b 

where M a b(x,y) = [(j) a (x), Xb(y)\c is the matrix defined in Eqs. (174]) and d75j) . We recall that this 
matrix has the form 


( 0 

-1 

0 

o N 

1 

u 

5 

T 

0 

V 

A 

B 

v° 

w 

c 

v) 


(97) 


where in the SU(ri) gauge field case, each entry in M is a 2 n x 2 n matrix. Using the block inversion 
method given in Eqs. (1A18I) and (1A19D . we find that M~ l is given by 


( E 1 

- 1 0 
TV + GW 0 
yUV + lW 0 


-(ST + TH) 
0 
T 
U 


-(.SQ + T1) ^ 
0 
G 

1 ) 


(98) 


where 


E =U-S(TV + GW) -T{HV + 1W) , 


(99) 


and where T, G , H, I are the elements of the block inversion of the matrix N of Eq. () 76 [) . 



( 100 ) 
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Substituting these into Eq. (1961) we find for the Dirac bracket a lengthy expression, which sim¬ 
plifies considerably after noting that [F(T),xi]c = [-F(T), —= 0 and [</>i, G'('h, d^)]c = 
flb)]c = 0, leaving the relatively simple formula 

[F(*),G($, * f )b =[F$), G ($, #)]c 

Xs]c (f [<h,G($, &)] c + G [04, G($, * f )]c) 

-[F$), X4]c (n [03, G($, * f )]c + I [04, G(#, * f )] C ) • 

( 101 ) 


We note that just as in the path integral derivation, only the matrix N enters, in this case through 
its inverse, rather than the full matrix of constraint brackets M. The final step is to evaluate the 
inverse block matrix elements J-, G, T~L, Z from the expressions for A, B , C, T>, again by using the 
block inversion formulas of Eqs. ([A18D and (IA19D . Let us define the Green’s function D _1 (x — y) 
by 

(i(L s ) 2 + a ■ LzxL s )T>- l (x-y) = 8*(x-y) , (102) 


and a second Green’s function Z{x — y) by 
Z(x - y)=A- BV~ l C 

= - 2 iga ■ B6 3 (x - y ) - 4 D s • -y)L$- *Dy 


(103) 


where we recall that for Tq = 0 gauge L = B + a x E. Then the needed inverse block matrices are 


G = - Z~ 1 BV~ 1 , 

H = ~V~ l CZ~ l , 

Z =V ~ 1 + V^CZ-'BV- 1 


(104) 

We wish now to apply the Dirac bracket formula to the cases (i) E('I') = T and 4H) = \k', 
and (ii) F(Jff) = T and G(T,T^) = H, with H the constrained Hamiltonian of Eq. (195|) . The 
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following canonical brackets are needed for this: 

[^(®),X3(i/)]c =2Dg6 3 {x-y) , 

['h(x),X 4 (y)]c ={iLg-a x Lg)5 3 (x-y) , 

[fa (x),&(y)]c =2 Dg5 3 (x-y) = -25 3 (x-y) < Dy , 

[0 4 (5'),'h T (y)]c = - ( iLg-Lg x a)d 3 (x- y) = S 3 (x - y)(iL$ - x a) , 

[fa (x),H] c =igB{x) • #(x) , 

[^4p), ff]c =^4 - 4 x a) x 4 ■ $(f) . 

(105) 


Additionally, for case (i) we need the canonical bracket 

[^i(^),^(y)]c =[^i(x),iPj(y) - e jk iP k (y)ai]c 

= - i(6ij + i€jnai)5 3 (x — y) = -iajaA 3 ^ - y) = -2 i(Sij - ^a i a j 'j5 3 (x - y) , 

(106) 


and for case (ii) we need the canonical bracket 

1 


[’k ?; (x), H] c = - [iDg x \P(x) - ax (Dg x \&(x))^ 


(107) 


Up to this point, we have not specialized L so as to make it easy to ascertain what the formulas 
become when g a Ug e fixi „ g is „ mi «ed (« in Q end [dj]). When L = 0, the matrix N degenerates 
to its upper left element A, which is a local function of x and so is algebraically invertible. For the 
Dirac bracket of IPj(x) with 'hj (y) we then find 

[^i(x),^t(^)] D =[^j(x),’p](y)]c - J z)[fa(z),^](y)]c 


= -2 i 
= — 2 i(x 


(fa - -a.a^^ix-y) - Dgi 6 15 Vj 

' Z / nrr • 


&ij 1 -)- 11 / _* Hj 

2 V ga-B 


ga ■ B(x) 

\y) , 


(108) 


where in the final line we have written j = ITj to relate to the abstract operator notation of Velo 
and Zwanziger 4], Multiplying the final line by i to convert the Dirac bracket to an anticommutator, 


and by a factor 1/2 reflecting our different field normalization, Eq. (11081) becomes the expression for 
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the anticommutator given in the zero mass limit of Eq. (4.12) of 
A, one can verify (as in Appendix C of jj]) that 


a. 


Using identities in Appendix 


((7 X Dx)i 


c _1 
0ij 




)*><*-a 

J ga ■ B(x) 


= 0 


(109) 


that is, the constraint x is explicitly projected to zero. 

Now setting L = B + a x E, so that the ui constraint implies To = Tq = 0, we find for the 
Dirac bracket of 'lq(x) with the constrained Hamiltonian 


* =[^(x),H] d = if)g x T(x) -ax (D s x T(x))] 
- / d 3 y^2Dg F(x - y)igB(y) ■ 'k(y) 

+Q(x - y)^(iLg- L$ X a) X Dg • $(y) 

+(iLx — ax Lg) H{x-y)igB(y) ■ 'h(y) 

+l(x - y)^(iL$-Lg x a) x Dg ■ $(y) | 


( 110 ) 


The first line of this equation is the unconstrained equation of motion in the form of Eq. m 
(when Aq = \ko = 0), while the remaining terms guarantee that 

#3 _dx _ d(a x D ■ T) _ -> d$_ 

dt dt dt dt ’ 

# 4 _X ■$ _ f d $_ 

dt dt ' dt 

( 111 ) 


where we have used the fact that we are assuming that D and L are time independent. This 
restriction can be avoided by treating the gauge fields as dynamical variables, taking into account 
their own constraint structure, and noting that the radiation gauge fixing constraint V • Py = 0, 
with the canonical momentum conjugate to A, has nonvanishing fermionic brackets with all 
Rarita-Schwinger constraints involving D = V + gA. This requires an extension of the Dirac 
bracket construction to take the new, Grassmann-odd, brackets into account, and the extended 
Dirac bracket structure will then obey Eq. (jllll) without requiring the assumption of a time 
independent A and L. 

In the next section, where we show positivity of the Dirac anticommutator, we will use covariant 
radiation gauge instead of 'I'o = 0 gauge. Setting L = D and putting everything together, we find 
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for the Dirac bracket of v f , *(x) with v I / ] (y), 

[^i{x),^]{y)\ D = - 2 i(Sij - ^0i(rj)6 3 (x - y) 

iT~{x yj 2D^iQ(x y jQy X ^ 

‘\- c 2(j l Drji a x Dx^ihL(x y)*Dyj (J'-^x cr x Dy^){nD y ^Ty x ; 

( 112 ) 

which gives the generalization of Eq. (11081) to the case when a covariant gauge fixing constraint is 
imposed. 


IX. QUANTIZATION OF THE ANTICOMMUTATOR DERIVED FROM THE DIRAC 
BRACKET AND POSITIVITY IN COVARIANT RADIATION GAUGE 


Given the Dirac bracket, the next step is to quantize, by multiplying all Dirac brackets by i and 
then reinterpreting them as anticommutators or commutators of operators. In the case considered 
here, this can be done in a constructive way, as follows. First let us replace the set of 2 n component 
column vector constraints <j) a and 2 n component row vector constraints Xa by the set of 4n scalars 
given by their individual matrix elements. Moreover, since the Xa are the adjoints of the 4> a , we can 
take linear combinations to make all of these scalars self-adjoint. Labeling the set of self-adjoint 
scalar constraints by <b a , the Dirac bracket construction for the bracket of F with G reads 


[F, G] D =[F, G\ c - E E^’ cT-^h, G\c 

a b 

Tab =[$a, $b]c , 


(113) 


with the matrix T real. 

We now observe that since the <J> a are all linear in the scalar components of T and \M, if we 
make the replacement i\ , \c —)► { , }c, with { , } the anticommutatior, and replace all Grassmann 
variables T and with operator variables having the standard canonical anticommutators, then 
since there is no other operator structure the same real matrix T a b will be obtained. Moreover, if 
F and G are both linear in the scalar components of T and \M, the Grassmann bracket i[F, G\c 
formed from scalar components of F and G will agree with the canonical anticommutator i{F, G}c 
formed from the corresponding operator scalar components, and will be a c-number. Thus, for 
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linear F and G we can define a “Dirac anticommutator” {F, G}d by 

{F, G}d ={F, G}c - E *a}cT- b X {^ G}c 

a b 

T a b={®a,®b}c ■ 


(114) 


When one or both of F and G is bilinear, the Grassmann bracket i[F, G\c formed from the scalar 
components of F and G will agree with the canonical commutator formed from the corresponding 
operator scalar components, and we can define a “Dirac commutator” by a formula analogous 
to Eq. (11141) in which each anticommutator with at least one bilinear argument is replaced by a 
commutator. In this way we get a mapping of classical brackets into quantum anticommutators and 
commutators, that inherits the algebraic properties of the Dirac bracket, including the chain rule, 
with the Jacobi identities for odd and even Grassmann variables mapping to the corresponding 
anticommutator and commutator Jacobi identities. 

To complete this correspondence, we must show that the Dirac anticommutator of and 
Ty 01 (with a = 1,2, j3 = 1,2 the spin indices, u = l,...,n, v = 1 ,...,n the internal symmetry 
indices, and i = 1,2, 3, j = 1, 2,3 the spatial vector indices) has the expected positivity properties 
of an operator anticommutator, by showing that for an arbitrary set of complex functions Af u (x ), 
we have 


d 3 xd 3 yAr(S)Af v (y){^r^)^V V (y)}D > 0 


(115) 


We demonstrate this in several steps, in covariant radiation gauge. First we examine the conditions 
for positivity of the canonical anticommutator and Poisson bracket, 

J d 3 xd 3 yAr(x)A*/ v (y){^r(x), v (y)} c = f d 3 xd 3 yAr(S)Af v (y)i[^f u (x), v (y)] c . 

(116) 

From '1'^ ;r = iPj V — ejkiP^crf 3 , we find that 


[*r0?), *Y v mc = - i{SijS al3 + ie jlk af)5 uv 5 3 (x-y) 


1 


= - i(aj(Ti) a,3 5 uv 5 3 (x -y) = - 2 *(<% - -<Jiaj) al3 6 uv 6 3 (x - y) 


(117) 


Multiplying by i/2, and writing Af u = Rf u +i = 1,2,3, a = 


1,2 , u = 1,..., n, with R and 


I real, the right hand side of Eq. (11161) evaluates to (we suppress the internal symmetry index u 
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from here on, so (Rf) 2 means Ylu=i(R? u ) 2 > etc.) 

EE((^ 2 +(^n 2 )-^((^2-4 1 +4 2 ) 2 +(^+4-^i) 2 +(^2 2 +4 2 +4 1 ) 2 +(^-4 2 +^3) 2 ) • 

j=l a=l 

(118) 

If all three components Af, i = 1,3 are present, the expression in Eq. (11181) is not positive 
semidefinite. But when only two of the three components are present, as a result of application of 
a constraint, then each of the four squared terms on the right hand side of Eq. (|118l) contains only 
two terms, and so the expression in Eq. (1118D is positive semidefinite by virtue of the inequality 


X 2 + Y 2 - ^(X ± Y) 2 = \{X^ Y) 2 > 0 


(119) 

Another way of seeing this, noted by both Velo and Zwanziger 4j and Allcock and Hall Q] , is 
that because X)jLi a i a i = 3, the expression W t] = 5ij — is not a projector. But when one 

component of a, say < 73 , is replaced by 0 , so that one has a i a i = Yli =1 = then 

1 1 2 ' 1 

y, WuWij = Sij - 2 + -ai y (JiCTiCTj = % - = Wij , (120) 

1 1=1 

and Wij is a projector and hence is positive semidefinite. So we anticipate that proving positivity 
will require projection of Eq. (11171) into a subspace obeying at least one constraint on 

The next step is to use the property that the Dirac bracket of linear quantities F and G reduces 
to the canonical bracket of their projections into the subspace obeying the constraints, when (as 
is the case here) all constraints are second class, that is they all appear in the Dirac bracket [181 ]. 
Referring to Eq. (I113D . let us define 

a b 

g=g-YY.I°-^cT^ , 


( 121 ) 


so that 


[F, <3? c ]c =[F, <& c ]c - y y [F, ^alcT-^b, <E 

a b 

=[E,$ c ] c -yy[F,$ 0 ] c r- i T bl 

a b 

=[F,d. c ] c -y[F,d. a ] c( 5 ac = 0 , 


c\C 


( 122 ) 
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and similarly for G. As a result of this relation, which holds when the canonical brackets are 
simply numbers (as in the case here where 4> c and F, G are linear), together with symmetry of the 
canonical bracket [G, $ c ]c = [3?o G\c, we see that 

[F,G} d = [F,G}c • (123) 

These properties of Eqs. (112111 (112311 carry over when we replace Grassmann numbers with opera¬ 
tors, and classical brackets with anticommutators, since in the linear case all anticommutators of 
linear quantities are c-numbers that commute with the operators, and since the anticommutator 
is symmetric. Thus we have 

. (124) 

To further study the properties of 40 (x) and T|(y) (with spinor indices suppressed), let us now 
return to our original labeling of the constraints by 0 a and x.a as in Eq. (I101D . so that we have in 
the Dirac bracket formalism 

Ti(x) = 4q(x) - ^^[^(z),x a ]cM~V b , (125) 

a b 

and a similar equation (with the roles of 0 a and \a interchanged) for 'kt(y), with a, b summed from 
3 to 4. We now note two important properties of this equation. The first is that it is invariant 
under replacement of the constraints Xa by any linear combination x' a = XbKba , with the matrix K 
nonsingular, since the factors K and K~ l cancel between x! a an d Tf( r 1 . (More generally, the Dirac 
bracket is invariant under replacement of the constraints by any nonsingular linear combination 
of the constraints, reflecting the fact that the Dirac bracket is a projector onto the subspace 
obeying the constraints, and this subspace is invariant under replacement of the constraints by 
any nonsingular linear combination of the constraints.) The second is that if we act on 40 (x) with 
either or (<r x Dg)i, we get zero. For example, recalling that in covariant radiation gauge 
D Si * i{x ) = 04 (x), we have (with spatial variable labels x suppressed) 

Di^i = 04 - ^2 XalcM^-fo = 04 - ^2 ^2 M4aM ab 0 & = ^~^2 ^ Ab ^ b = 0 ’ ( 126 ) 

a b a b b 

and similarly for (er x D^)i, with 0 4 replaced by 03 . 

Let us now write 'Lj(x) as a projector Rij(x,y) acting on 'Lj(y), giving after an integration by 
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parts on y, 

Vi(x) = J d 3 yR ij (x,y)'if j (y) , 

Rij(x, y) =5ij5 3 (x -y) + EE/ d 3 z[^i{x),Xa{z)} c Mj(z,y)^ bj {y) , 

a b J 

(127) 

with 

V 3 j(y) = (o' X , V4j(y) = *Dyj • (128) 

By virtue of Eq. (11261) and its analog for a x D, we have 

DgiRij^x^y) =0 , 

(a x D s )iRij(x,y) = 0 , 

(129) 

and by the reasoning of Eqs. (1401) (|42|) we also have (assuming a ■ D is invertible) 

<TiRij(x,y) = 0 . (130) 

Next let us focus on the bracket [4h(x), Xa{z)\c appearing as the first factor inside the sum. Setting 
L = D in Eq. (11051) we have 

[^(£),X 3 {z)\c =2D s 5 3 (x - z) , 

[$(x),X4(z)]c =(iDg ~ a x Dg)8 3 (x - z) . 

(131) 

Using the invariance of 4q:, or equivalently of Rij, under replacement of X'i- X 4 by an y nondegenerate 
linear combination of X 3 , X 4 > let us choose the new combinations so that 

[^0?),X3(u)]c =(<? x Dg)S 3 (x — z) = ff 3 (x)5 3 (x- z) , 

[#(x), X4{z)\c =Dx^{x - y) = fj 4 (x)8 3 {x - z) . 

(132) 

Substituting this into Eq. (11271) . we get the symmetric expression 

Rij{x,y) = 8ij8 3 {x-y)+ / d 3 zfj a (x)M^ (x, y) V b j (y) 

a b ^ 


( 133 ) 
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By virtue of this symmetry, the projector R.y is annihilated by the constraints Uyj and (<? x D v )j 
acting from the right, which in turn implies that in addition to Eq. (11301) we also have 


Rij(x,y)aj = 0 


(134) 


An explicit construction of Rij(x,y) and verification of Eqs. (11301) and (11341) is given in Appendix 

C. 

Returning now to Eqs. (11151) and (11241) . writing T" and in terms of projectors acting on 
V P“ and \Ht we have (using = < 7 *f e , and continuing to suppress internal symmetry indices u, v, 
which are contracted in the same pattern as the spatial vector and spin indices) 


J d 3 x J d 3 yAf (x)A* 0 {y)(x ), T] lj (y)} D 
Jd 3 xJ d 3 yA?(x)A*\y){*?(S),&/(y)} c 

J d 3 x J d 3 yAf(x)A*d(y) J d 3 z J d 3 w R^(x, ^J(w)} c R*^(y, w) 

J d 3 x J d 3 yA?(x)A*d(y) J d 3 z J d 3 w R^[x, z)2 (di m S lS - 5 3 (z - w)R*^(y,w) 

2 / d?>Z / z) J d 3 yA 13 (y)R^(y, z) 

(135) 


which is positive semidefinite. 

We conclude that the anticommutator of T with 4U is positive semidefinite in covariant radiation 
gauge. The symmetry of the ^3 4 and X 3,4 constraints in this gauge is essential to reaching this 
conclusion; if gauge fixing were omitted, or if another gauge were chosen, this symmetry would not 
be present and we could not deduce positivity in a similar fashion. 


X. LORENTZ COVARIANCE OF COVARIANT RADIATION AND T 0 = 0 GAUGE 
AND LORENTZ INVARIANCE OF THE DIRAC BRACKET 

We study next the behavior of covariant radiation gauge and the Dirac bracket under Lorentz 
boosts. The Rarita-Schwinger field and its left-handed chiral projection T" both have a four- 
vector index y and a spinor index a. Under an infinitesimal Lorentz transformation, the trans¬ 
formations acting on these two types of indices are additive, and so can be considered separately. 
The spinor indices are transformed as in the usual spin ^ Dirac equation by a matrix constructed 
from the Dirac gamma matrices, which commutes with D^. Hence the spinor index transformation 
leaves the covariant radiation gauge condition D ■ T invariant. 
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This leaves the transformation on the vector index to be considered, and this is a direct analog of 
the Lorentz transformation of radiation gauge in quantum electrodynamics [3] . Since the radiation 
gauge condition is invariant under spatial rotations, we only have to consider a Lorentz boost, 


x ~^X = x + vt , 
x° =t —> t' = t + v ■ x . 

(136) 


Under this boost, the field T transforms as 

T —> T 7 = T + FT 0 


(137) 


For an observer in the boosted frame, covariant radiation gauge would be Dgi ■ T 7 = 0, with 
Dtfi = Vg> + gA', where A' = A + O(u). Applying this to T 7 (x 7 ,t 7 ) and using the covariant 
radiation gauge condition in the initial frame, we get 

Dgi ■ T 7 = VjT,j(x,t ) , (138) 

with S j(x,t) a local polynomial in 'L, To and the gauge fields, where we have dropped primes on 
the right hand side since there is an explicit factor of v. So in the boosted frame T 7 does not obey 
the covariant radiation gauge condition, but this can be restored by making a gauge transformation 

T 7 -> — D(D 2 )~ 1 VjT, j (x, t ) . (139) 


Hence the covariant radiation gauge condition is Lorentz boost covariant, although not Lorentz 
boost invariant. Analogous arguments apply to To = 0 gauge. Under the Lorentz boost of Eq. 
(11361) . T 0 —> T 0 — v- T, so T 0 = 0 gauge is not preserved, but it can be restored by a compensating 
infinitesimal gauge transformation. 

Referring now to Eq. (1C10I) . we note that the covariant radiation gauge Dirac bracket and 
the anticommutation relations are invariant under infinitesimal Rarita-Schwinger gauge transfor¬ 
mations, such as that of Eq. (11391) . up to a remainder that is quadratic in the gauge parameter. 
Hence the covariant radiation gauge Dirac bracket and the anticommutation relations following 
from it are Lorentz invariant, since a finite Lorentz transformation can be built up from a series of 
infinitesimal ones. 
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XI. COMPATIBILITY OF LORENTZ COVARIANCE AND MODE COUNTING IN 
ON-SHELL RARITA-SCHWINGER PARTICLE-PHOTON SCATTERING 

We address finally the question Q of whether leading order on-shell scattering of Rarita- 
Schwinger fields from an external electromagnetic field has the requisite relativistic covariance, 
while preserving the correct counting of massless spin 4 propagation modes. The operator effective 
action for this scattering process can be read off from the interaction term in Eq. CD, 

£ e s(iPn,A v ) =^g ip ll (x)ie firiup 'y 5 'y 71 A u (x)'ip p (x) , 

(140) 

where we have suppressed spinor indices as in the text from Eq. © onwards. Since we are taking 
the external field A v as Abelian, the covariant derivatives in the equations of motion and constraints 
are given by 

D u = d u + gA u , *D l/ = *d u - gA v . (141) 

To keep the analysis simple, we shall assume that A v [x) is of short range, and vanishes for \x\ > 
R for some radius R. This effective action, the equations of motion of Eqs. © and ©, and 
the primary and secondary constraints following from them, given in Eqs. © and ©, are all 
relativistically covariant. Hence to do a covariant calculation, it is only necessary to choose a 
relativistically covariant gauge fixing, which for convenience we take as the gauge covariant Lorentz 
gauge condition 

D p % = 0 , (142) 

which is attainable from a generic gauge by the gauge transformation of Eq. ©, provided that 
D p D p is invertible. 

Taking the matrix element of Eq. (j!40j) between an incoming Rarita-Schwinger state of four- 
momentum p. and an outgoing Rarita-Schwinger state of four momentum p' , we get the corre¬ 
sponding scattering amplitude 

A= \ i9 J d 4 x ^v(p'’ x ^ eMUP ^v A v( x )^p(Pi x ) > ( 143 ) 

where and are now wave functions, rather than operators, that obey the Rarita-Schwinger 
equations of motion in the presence of the external held A v . 
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We now introduce source currents for the gauge potential A u and the Rarita-Schwinger wave 
functions if> p and and study their conservation properties. The source current to which the 
gauge potential A v couples is defined by writing the scattering amplitude as 

A =— ig j d 4 xA l/ (x)J u (x) , 

J v (x) =V’ /i (p',®)e /ir ' i/p 757 n^ P ip,x) . 

(144) 

The source current for the Rarita-Schwinger field 'ifj.Jp', x) is defined by writing the scattering 
amplitude as 

1 


A= 2 i S J '^x^iAp, x)J ll (p, x) , 
J' l (p,x) =e mup 'y 5 'y v A u (x)ifj p (p,x) . 


(145) 


Finally, the source current for the Rarita-Schwinger field ip p (p, x) is defined by writing the scattering 
amplitude as 

■A=\ig J d 4 xj p (p' ,x)ipp(p,x) , 


J AP',x) =ip fl {p',x)e priup 'y 5 -f v A u (x) . 


(146) 


We now show that the three currents that we have just defined are conserved. For the source 
current J u for the gauge potential, we have 

duJ v =Vv(lA x)5 v e mvp ''%'y rt '4>p(p,x) 

x)e pr]Vp ^ v D v \l) p {p, x) 

=0 , 

(147) 

where the first and second terms on the right vanish by the Rarita-Schwinger equations for 'ipAp 1 , x) 
and ip p (p,x) respectively. For the source current J p pp,x) for the spinor ipAp^x) , we have 

DpJ p {p, x) =e W(/p 757 r, {dpA v {x))ij}p[p, x) 

+e M??!/p 757 v A u (x)Dp'if>p(p, x) 

=0 . 


( 148 ) 



39 


The second term on the right vanishes by the Rarita-Schwinger equation for ijj p (p,x), while the 
first term on the right can be rewritten as 

^ vup 'y 5 -f v F^(x)ilj p (p,x) ( 149 ) 

and vanishes by the secondary constraint of Eq. JH]). Finally, for the source current J P (p\x) for 
the spinor ij) p (p,x), we have 

J P (p, x)J5 p =^ p {p\ x)e^ up 'y 5 'y v (' d p A v (x )) 

+^ p ip\ X )*Dpt pwp lslr 1 Av{x) 

=0 . 

( 150 ) 

Again, the second term on the right vanishes by the Rarita-Schwinger equation, while the first 
term on the right vanishes by the secondary constraint of Eq. Q . 

Consider now the following three gauge transformations, 


A v (x) — >A v (x) + d u A , 

^p(p,x) -+i> p (jp,x) + D p a , 
x) -^(p', x) + pEp , 

(151) 

with a and f3 independent spinorial gauge parameters. From Eqs. (|144jl - (|146jl . together with Eqs. 
(|147jl - (jl50jl . we find that these transformations all leave the amplitude A invariant, 

5\A =—ig j d i x(d u A)J u (x) = —-ig J d 4 xAd u J u (x) = 0 , 

5 a A =ig J d 4 xJ p (p,x)D p a =-^ig J d 4 xj p (p, x)J5 p a = 0 , 
deA =Ug J d A x!35 p J p (jp,x) = ~ig j d‘ i x/3D p J p (p,x) = 0 . 

(152) 


We next must specify more precisely the structure of the spinor wave functions entering the 
formula for A. Since the gauge field A u is assumed to vanish in the external region \x\ > R, the 
Rarita-Schwinger wave functions obey free field equations in this region. So for \x\ » R they can 
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be taken asymptotically as plane waves at t ±oo, 

ipp(p',x) ~ Ufl (p')e ip '' x , t->+ oo 
ipp(jP,x ) ~ Up {p)e ip ' x , t ->• -oo 


(153) 


With these boundary conditions, the formula for the amplitude takes the final form 

A = ^ig J ] (j)\x)e mup ^ v A v (x)^{p,x) . (154) 

The out state (-) and in state (+) boundary conditions used here are analogs of the boundary 


conditions used in the distorted wave Born approximation 


2l|, which the construction of Eq. (11541) 


resembles. Equation (11540 then gives an approximation to the matrix element for Rarita-Schwinger 
scattering by the gauge potential. 

Finally, using the gauge invariances of A, we can now address the question of counting degrees 
of freedom. We have imposed one gauge fixing constraint on the Rarita-Schwinger wave func¬ 
tions, given in operator form by Eq. (1142 1) . which implies that the wave functions obey analogous 
conditions of vanishing gauge covariant divergence, 


i’ Jp',x)l5 p = D p il>p{ Pl x) = 0 


(155) 


To see what further conditions can be imposed, we must examine consistency with the Rarita- 
Schwinger equation for the operator ipp. This is most easily done using the alternative form of the 
Rarita-Schwinger equation given in Eq. (1A6I) . 

7 p {D v ip p - Dpip v ) = 0 . (156) 

Applying D u from the left to this equation, we get 

D v D vl p i\> p - 7 p [D v , D p ]il> v - YDpD v ipv = 0 , (157) 

which using [D U ,D P ] = gF vp becomes 

D v D vl p ^ p - 57 P F vp ^ ~ l P D p D^ u = 0 . (158) 

When = 0, this equation simplifies to 


D v D vl p i\) p = g lp F vp ^ v 


(159) 


In regions where the gauge held is nonvanishing, the right hand side of this equation is nonzero, 
and so we cannot impose the additional constraint 7 p ijjp = 0. This is not surprising, since where the 
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fields are nonzero there is already a second covariant constraint given by the secondary constraint 
of Eq. (J 6 ]), which in terms of the dual field F up = i e vpa ^F a ^ takes the form 


0 = 7 p F^ u , 


(160) 


which differs in form from the right hand side of Eq. (11591) when the fields are not self-dual or 
anti-self-dual. However, in the external region where the gauge field vanishes, Eq. (I159D becomes 


d v d vl p ^ P = 0 


(161) 


and we can then impose 'y p ip p = 0 by making a gauge transformation in the external region 
that preserves the original constraint d p ip p = 0. This is possible, since the condition for the 
transformation ip p —>• t/jp + d p a to preserve d p p' p = 0 is d p d p a = 0 , which is compatible with the 
form of Eq. (I161D . This establishes that in the external region we can impose the two Lorentz 
covariant conditions d pr p) p = 0 and 'y p ipp = 0, which together with the primary constraint of Eq. 
(|5|) suffices to give the correct degree of freedom counting for the incoming and outgoing Rarita- 


Schwinger wave functions 


22], Alternatively, if we are not concerned to maintain manifest Lorentz 


covariance, we can make a gauge transformation in the external region to the gauge i/jq = V • ip = 0 


used in 


231 ] to enumerate Rarita-Schwinger degrees of freedom. 


Note that if one were to construct the Born approximation amplitude, in which the Rarita- 
Schwinger wave functions in the presence of the gauge field are replaced by plane waves in the 
interior region where the potential is nonzero, the arguments given above for compatibility of 
Lorentz covariance with degree of freedom counting would fail. The reason for this is that the 
spinor source currents would then no longer be conserved, even to zeroth order in the gauge 
coupling g, because the free particle plane wave solutions do not obey the secondary constraint 
of Eq. ([ 6 ]) . The non-existence of a satisfactory Born approximation for Rarita-Schwinger photon 
scatterin g ag rees with restrictions on massless particle three point functions obtained by “on-shell” 


methods 


24]. It is also related to the fact that the extra gauge terms d u £F in Eq. (5.27) of 23] for 


the free Rarita-Schwinger propagator are an obstacle to setting up a Lippmann-Schwinger equation 
for Rarita-Schwinger photon scattering, which could be expanded in a perturbation series in g. To 
establish compatibility, we have had to use an analog of the distorted wave Born approximation 
21 ], in which the leading approximation to the amplitude is constructed using interacting rather 
than free fermion wave functions and does not have a gauge-invariant small coupling, g —> 0 limit. 

The results just obtained will play a role in the analysis of closed Rarita-Schwinger fermion 
loops with gauge field vertices. We conjecture that when fermion propagators corresponding to 
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the interacting wave functions are used in the loops, the usual vertex conservation properties 
including the known chiral anomalies will hold, while permitting changes of fermionic gauge of the 
Rarita-Schwinger fermions in the loops. We plan to pursue this as a further investigation. 


XII. DISCUSSION 


To conclude, we see that unlike the massive case, the massless gauged Rarita-Schwinger equation 
leads to a consistent theory at both the classical and quantized levels, and these statements are 
already foreshadowed in the equations of [4] when taken to the zero mass limit. Thus, non- 
Abelian gauging of Rarita-Schwinger fields can be contemplated as part of the anomaly cancelation 
mechanism in constructing grand unified models. 

Our analysis invites a number of extensions: 


1. We have derived the path integral formula and the constrained Hamilton equations of motion 
in To = 0 gauge, and shown positivity of the Dirac bracket in covariant radiation gauge. This 
is reminiscent of what happens in quantizing non-Abelian gauge theories, where unitarity is 
manifest in one gauge, and renormalizability in another. In the gauge theory case, one has 
an apparatus for transforming from one gauge to another, but that remains to be worked out 
for fermionic gauge transformations of the Rarita-Schwinger fields. The analysis of Sec. 11, 
where we computed an approximation to Rarita-Schwinger photon scattering and found full 
fermionic gauge invariance, is evidence that at least on shell, where the equations of motion 
are obeyed, different gauges are equivalent. 

2. The analysis of Sec. 11 shows that to get a nonvanishing approximation to Rarita-Schwinger 

photon scattering, one has to go beyond the Born approximation to scattering, by using in 
and out solutions of the scattering equations in the presence of the photon perturbation. (We 
have assumed the existence of these interacting solutions of the Rarita-Schwinger equations, 
but this remains to be studied.) This is a consequence of the fact that the secondary 
constraint u = 0, although involving the gauge fields, is zeroth order in the gauge held 
coupling g because it arose from application of to the Rarita-Schwinger equation 

of motion. Thus a quantum theory of Rarita-Schwinger fields in a gauge held background 
necessarily has a non-perturbative aspect. 

3. A possible exception to the non-perturbative behavior of fermion lines is when the E and B 



43 


gauge fields are random, since if Eq. (1501) is replaced by an average, denoted by AV, 

<*»>“ - (^) AV ' ■ < 162 > 

it becomes 


(To)av — <? • Wav , 


(163) 


which is compatible with (Tq)av = ■ (T)av = 0, the customary free Rarita-Schwinger 


constraints employed in 


23|. This heuristic observation suggests that Rarita-Schwinger fields 


coupled to quantized gauge fields with zero background gauge field may have a perturbative 
g —> 0 limit. 


4. We have included in our analysis only fermionic constraints. However, the action exponent 
in Eq. (1511) has the form (in Aq = 0 gauge) 


1 

2 


crx 


— \E^ • (B + <7 x E)(a ■ B) V-DxT 


1 


-W • a x D(a ■ B) 1 (B + a x E) • T + T+ • D x T - x a ■ 


■ j d 4 x{d 0 $ ■ 


P — bosonic constraints) — / dtH , 


(164) 


with H the Hamiltonian obtained from the stress-energy tensor. We have not included the 
bosonic constraints in our analysis, and we note that they will have fermionic brackets with 
the fermionic constraints (j> a , Xa■ (For a discussion of bosonic versus fermionic constraints, 


see 


25 |.) 


5. In quantizing, we assumed that the gauge fields A are time independent, so that d/dt and D 
commute. As already noted, this assumption can be dropped if the gauge fields are treated 
as dynamical variables, leading to an extension, yet to be analyzed, of the bracket structure, 
again involving fermionic brackets. 


6. In showing in the Abelian case that there is no superluminal propagation, the inversion of 
a ■ B to get To only required ( B ) 2 ^ 0. In the non-Abelian case, where B is itself a matrix, 
the conditions for invertibility are nontrivial and have yet to be analyzed. 


7. In demonstrating positivity of the anticommutator, we made essential use of the condition 
(T ■ T = 0. Deriving this from the covariant radiation gauge condition D ■ T = 0 assumed the 
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invertibility of 3-D, and attainability of covariant radiation gauge assumed the invertibility 
of (-D) 2 . The conditions for invertibility of these two operators remain to be studied. (The 


open space index theorems of Callias 


26j ] and Weinberg [26y involve 3 ■ D + iq i>, with cj) a 


scalar field, and so do not give information about the invertibility of 3 ■ D.) 


8 . The fact that when gauge fields are present, the gauge conditions To = 0 and D ■ T = 0 are 
not equivalent, has consequence for the study of condensates that give rise to a mass for the 
initially massless Rarita-Schwinger field. The natural way to form a condensate involving the 
spin | field r/v an d a spin | field A is through a structure of the form If To and 3 ■ T 

were simultaneously zero, this would translate back in the covariant formalism to 7 m V ; m = 0 ) 
forbidding such condensate formation. Thus, the possibility of dynamical mass generation 
for the spin | field depends on the inequivalence that we have demonstrated of the gauge 
conditions To = 0 and D • T = 0 in nonzero gauge field backgrounds. An alternative way 
to see that the gauge condition = 0 cannot be imposed when fields are present is to 

note that its left chiral reduction is To = 3 • T, which in general conflicts with the secondary 
constraint 3 ■ B T 0 = (B + 3 x E) • T. 
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Appendix A: Notational conventions and useful identities 


We follow in general the notational conventions of the book Supergravity by Freedman and 
Van Proeyen 23]. The metric r ] fll/ is (—, +, +, +) and the Dirac gamma matrices 7 ^ , 7 ^ obey the 


Clifford algebra 


7m T' + 7^ 7m = 2 7m^ 


(Al) 
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They are given in terms of Pauli matrices 07 by 


7o = -7 = 


Ifj = l 3 = 


75 = *70717273 = 


0 -1 
1 0 


0 a 




°3 0 


1 0 
0 -1 


We also note that 


, _ ,0123 _ 1 

€0123 — — € — 1 , 


the left chiral projector Pi is given by 


-Pl =-(I + 75) , 


and the spinor ip is defined in terms of the adjoint spinor ip' by 

ip = ip'i'y 0 


(A2) 


(A3) 


(A4) 


(A5) 


As noted in |23(], the Rarita-Schwinger equation of motion can be written in a number of 
equivalent forms. When ordinary derivatives are replaced by gauge covariant derivatives, these are 
the vector-spinor equations 

e^ p lv D u iP p =0 , 
l wp D v 'P p =0 , 

7 p (D u ip p - D p ip v ) =0 , 

7 a D a {D a ip v - D u ip a ) =y([D p ,D rT \ip„ + [D u ,D p ]ip a + [D a , D u ]ip p ^j , 

(A 6 ) 

with only the fourth line, which is quadratic in the covariant derivative, involving more than just 
a substitution d v D u in the formulas of [23j. Using r y v y up = 2y p , these also imply the spinor 


equation 7 up D u ip p = 0. These formulas play a role in verifying stress-energy tensor conservation, 


as does the identity 


M 


0 — e' K<JIW {A T B\C rJ D p E y + A v B r C\D,jE p + A p B u C T D\E,j + A a B p C u D r E\ + AxB a C p D u E T 


(A7) 
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with A t , B\ , C a , DE u five arbitrary four vectors. This identity follows from 

0 = 8^e Xa ^ + 5 v T e aXatl + 5^e vaXa + 5"e^ aX + 5 x e a ^ a , (A8) 

which is easily verified by noting that A, <r, v must take distinct values from the set 0,1,2,3, 
and that r must be equal to one of these values. 

The fundamental identity for the Pauli matrices is 

^ah T itabc&c , (A9) 

with ei 23 = 1 and with the index c summed. We repeatedly use the following two identities that 
can be derived from Eq. (1A9I) . for a general three vector A that is proportional to a unit matrix 
in the spinor space and so commutes with 3, 

3 x (<r x A) = — 2 A + iff x A , 

(A x 3) x 3 = — 2 A + iA x 3 . 

(A10) 

Additional useful identities are 


3 x 3 =2i3 , 

33 - A =A — ia x A , 

3 ■ A 3 =A + ia x A , 

{3 x A) ■ 3 = — 2i3 ■ A , 

3 ■ (3 x A) =2i3 ■ A , 

&a&b = 2(<^afe 2 ®b&a) i 

B = iA — A x 3 = —(B x 3) 

(All) 

Gauge field covariant derivatives are 

Dp = d^, + gA^ , (A12) 

with the gauge potential A p = A^a and the gauge generators anti-self-adjoint, and with the 
components A') self-adjoint. The non-Abelian generators tA obey the compact Lie algebra 


[tA,tB\ = fABCtc ; 


(A13) 
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in the Abelian case we replace tA by — i- In writing field strengths E and B we pull out an 
additional factor of i to make them self-adjoint, so that we have the identities 


D x D = — igB , 

[D, Do] = - igE . 

(A14) 

We will also write a right-acting three-vector covariant derivative as + gA, and define a 

left-acting three-vector covariant derivative as jD = ^ — gA, so that we have the integration by 
parts formulas 

J dPxAJ^^B = — J d 3 xAiDgB , 

D s 5 3 (x - y) = - 5 3 {x - y)Dg . 

(A15) 


<— 


An analogous definition is used for the operators L and L which enter the gauge fixing condition. 

At the classical level, variables will be either Grassmann even or odd. Irrespective of the 
Grassmann parity of monomials A and B, the adjoint operation is defined by 23] 


(AB)t = . 


(A16) 


For classical brackets, we follow the convention of Henneaux and Teitelboim 


28 ], 


,'dFdG dFdG . 
’ C ' dq L dpi dpi dq l ' 




d L F d L G d L F d L G 

+ 


dO a 9vr n 


dTT n d9 a 


(A17) 


with €f the Grassmann parity of F, with d L a Grassmann derivative acting from the left, and 
with q i , pi ( 9 a , ir a ) canonical coordinates and momenta of even (odd) Grassmann parity. Using 
the classical bracket, the Dirac bracket is constructed from the constraints as in Eq. (|96D of the 
text. To make the transition to quantum theory, the quantum commutator (anticommutator) is 
defined to be ih times the corresponding Dirac bracket (with h = 1 in our notation). Classical 
canonical brackets are always denoted, as above, by a subscript C, with a subscript D used for the 
corresponding Dirac bracket. We use the standard notations [A, B] = AB — BA for the commutator 
and {A, B} = AB + BA for the anticommutator. 
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To calculate the Dirac bracket, we use block inversion of a matrix. Let 

. Ai A 2 

M = 


M~ l = 



(A18) 


with Ai,...,A 4 themselves matrices. Then when A 4 is non-singular, the blocks Bi,...,B 4 of M 1 
are given by 


A 

=A\ 

— A 2 A 4 3 Az 

Bi 

=A~ 

-1 

5 

b 2 

= - 

A 1 A 2 A 4 1 

Bs 

= - 

a^a 3 A- 1 

B 4 

=*k 

_1 + a 4 3 a 3 a 


(A19) 


Even though the blocks are noncommutative, Eqs. (1A18I) and (1A19D give an inverse that obeys 
M~ l M = MM' 1 = 1 . 

When the constraints (j) a and \a are combined into an 8 element set of constraints K a = 
0m Ka+4 = Xa, cl = 1,..., 4 then the bracket matrix S a b(x,y) = [K a (x), Hb(y)]c can be expressed in 
terms of the matrix M a h ( x , y ) of Eq. d74|) as 


S(x,y) = 


0 M(x, y) 
M T (y,x) 0 


(A20) 


where Mj b (x,y) = Al^ a {x,y) is the matrix transpose. Defining the inverse M 1 (x,y) that obeys 
f d 3 zAf- 1 (x, z)M(z, y) = f d 3 zAI(x, z)M~ 3 (z, y) = 5 3 (x — y), it is easy to verify that 

0 M T ~ 1 (y,x) 

M -1 (x, y) 0 


S (x, y) = 


(A21) 


Appendix B: Analysis of the Rarita-Schwinger field in an external Abelian gauge field, 
without using the covariant radiation gauge constraint: propagation of the longitudinal mode 


We continue here the analysis begun in Sec. 5, but now without using the covariant radiation 
gauge constraint D ■ T = 0, to analyze propagation of the longitudinal mode that is eliminated by 
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gauge fixing. We must now solve for Cg"*" starting from Eq. (I64|) with Ci 2 = 0, so the third line 
of Eq. (jMj) simplifies to 

0=(B) 2 ndZ’±-K(Q 3 C 3 )^ , 

Q3 =B\E2 — B2E1 + B^cr ■ (B + iE ) — iB ■ Ea 3 

(Bl) 


Writing this as 

M / u u t/ 12 \ /elf 

V 0 / \ t/211/22 ) ^ ci j 

we find for the matrix elements 


(B2) 


Uu =(Bj 2 n - K[B x E 2 - B 2 E 1 - i(BiE 1 + B2E2) + B f] , 

U22 ~ K[B ] E 2 - B 2 E 1 + i{B x E 1 + B 2 E 2 ) - B |] , 

U 12 = — Af? 3 [i?i + iEi — i(f ?2 + iE/2)] , 

U21 = — Ai?3[Ai + iE\ + i(i?2 + i-E/2)] 

(B 3 ) 

The equation 0 = det(U) = U11U22 — U12U21 reduces, after dividing by an overall factor of (A) 2 , 
to 

0 = {B) 2 ^ 2 - 2 VtK(BiE 2 -B 2 E l ) + K 2 {El +El-Bl) , (B 4 ) 

with the solution 

n _a±t 1 /2 

K ~ 0)2 ’ 

X =B\E 2 — B2E1 , 

Y =(BiE 2 - B 2 E 1 ) 2 - ( B) 2 (E 2 + E 2 - B 2 ) . 

(B5) 

The analysis of the solutions of Eqs. ( 1 B 4 I) and ( 1 B 5 |) divides into two cases, according to whether 
the roots of Eq. (lB 5 j) are both real, or both complex. The roots are both complex if 

{B l E 2 - B 2 E l f < (B) 2 (E 2 + E 2 - B 2 ) , (B 6 ) 
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which can be rearranged algebraically to the form 

[(B) 2 - (E 2 + E 2 )]B 2 < (.B 2 + B 2 )(E 2 + E 2 ) cos 2 0 , (B7) 

where we have written 

B X E 2 - B 2 E 1 =(Bl + Bl) 1/2 (El + Elf' 2 sin0 , 

B l E l + B 2 E 2 =(Bl + B\) x I 2 (E\ + El) 1 ' 2 cos </> . 

(B8) 

Since the right hand side of Eq. (1B7I) is non-negative, when the left hand side is negative the 
inequality is satisfied, and both roots are complex. Hence a necessary (but not sufficient) condition 
for both roots to be real is 

(B) 2 - (El + E 2 ) > 0 . (B9) 

1. The hyperbolic case: both roots real 

When both roots are real, Eq. (jBlI) describes the hyperbolic case of propagating waves. Intro¬ 
ducing the velocity V = Sl/K, Eq. ()B4j) can be written as 

0 = (B) 2 V 2 - 2V(B 1 E 2 - B 2 Ei) + E 2 + E 2 -B 2 , (BIO) 

which can be rearranged algebraically to the form 

[(Bl + Bl) 1 / 2 -(El + El) 1 /' 2 ) 2 + (B) 2 (V 2 - 1) = 2(Bl+ Bl) 1 / 2 (E'l +El) 1 / 2 (V sin <p-l) . (Bll) 

Let us now assume that V 2 > 1, and show that this leads to a contradiction. When V 2 > 1, 
the left hand side of Eq. (IB 111) is nonnegative, which implies that V sin <p on the right must be 
nonnegative, and so can be replaced by its absolute value. Hence the right hand side of Eq. (IB111) 
obeys the inequality 

2(Bf+Bl) 1 / 2 (Ef+El) 1 / 2 (V sin0—1) = 2(B\+Bl) l/2 (El+El) 1 / 2 (\V sin (j)\-l) < 2(B) 2 (\V\-l) , 

(B12) 

where we have used Eq. (1B9I) . But the left hand side of Eq. (1B11D obeys the inequality 

[(B 2 + Bl) 1 ' 2 - (E 2 + El) 1 ' 2 ] 2 + (B) 2 (V 2 - 1) > (B) 2 (\V\ + 1)(\V\-1) > 2(H) 2 (|H|-1) , (B13) 

which is a contradiction, since a real number cannot be strictly less than itself. Hence we must 
have V 2 < 1, and there is no superluminal propagation. 
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2. The elliptic case: both roots complex 


When both roots are complex, Eq. m describes the elliptic case in which there are no 
propagating waves; when a propagating wave enters an elliptic region from a hyperbolic one it will 
be damped to zero amplitude. However, in the case of weak damping one can still define a wave 
velocity and ask what its magnitude is. When both roots are imaginary, Eq. (1B5I) takes the form 

n _x±i(-Y ) 1 / 2 
k~ (gy ’ 

X =B\E/2 — B 2 E 1 , 

-Y = - (BxE 2 - B 2 E 1 f + (B) 2 (E 2 + El - Bj) . 

(B14) 


Regarding H as real and the wave number K as complex, the effective propagation velocity has 
the magnitude 


IKd 


n 

~Er 


X 2 -Y E 2 +E 2 - B 2 
(B) 2 \X\ ~ \B 1 E 2 -B 2 E 1 \ 


The condition for weak damping is — Y « X 2 , which can be rewritten as 


(B15) 


(Bf(El + E 2 - B 2 ) « 2 {B x E 2 - B 2 E 1 ) 2 , (B16) 


and implies 


|Reff| « 


2\B\E2 — B 2 E 1 

W 2 


2| E 
\B\ 


(B17) 


Hence as long as 2\E\ is not much larger than \B\, which is required by the vacuum stability 
condition \E\ < \B\, the damped wave propagation velocity is subluminal. 


Appendix C: Construction of the projector Rij(x,y) 

Since there are only two constraints and two Xa constraints, we index them a = 1,2 rather 
than a = 3,4 as in the text, and use the invariance of Rij(x,y ) under changing the linear combi¬ 
nation of the Xa constraints. We start from the constraint set 

0! = a x D ■ # , xi = P -^5 , 

02 = • T , Xi = E ■ G ^ E . 


(Cl) 










For the bracket matrix 
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Mab(x,y) = [4> a {x),Xb{y)]c = 

we find the matrix elements 

A = — iga ■ B5 3 (x — y) , 

B =( 2 {D s f + ga ■ B)5 3 (x - y) = 5 3 (x - y)( 2(^) 2 + ga • B) 
C =(D g ) 2 8 3 (x - y) = 5 3 {x - y){5$) 2 , 

V =iga ■ B5 3 (x — y) . 



(C2) 


(C3) 


We write the inverse matrix M l (z,w) as 



which obeys 



(C4) 


(C5) 


In terms of the inverse matrix, the projector Rij(x,ui ) is given by (with internal symmetry indices 
suppressed) 


Rij(x,w) =5ij5 3 (x — w) 1 

~^~B)^iJ~{x rc)((T x “b B g ^Q{x w j A g j 

+(<T x D^iiiix - w)(a x ^DiAlj + ((? x D g )il(x - w)^5^j 


(C6) 


From this expression, we find 

D Si Rij(x,w) = Rijix^w^Drfj = (<r x Dz)iRij(x,w) = Rij(x,w)(a x *Drf)j = 0 . (C7) 

In verifying these, it is not necessary to evaluate the inverse matrix; instead, after contracting on 
the vector index i or j one expresses the resulting pre- or post- factor in terms of A, ...,T> and 
then uses the algebraic relations following from multiplying out the matrices in Eq. (jC5lh Finally, 
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contracting 


a ■ D s ai =(D S + ia x D s ) 
*Dw = (Dm — iff x 


a; a 


WJJ 5 


(C8) 


with Rij(x,w), we conclude that 


aiRij(x,y) = R ij (x,y)a j = 0 , (C9) 

when a ■ D is invertible. 

As a consequence of Eqs. (11271) and (1C7() . l'j(x) is invariant under the transformations 

T-^'F + .De , 

+ (j x De . 

(CIO) 


The first of these implies that the canonical anticommutation relations are invariant under in¬ 
finitesimal Rarita-Schwinger gauge transformations starting from covariant radiation gauge. 
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